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A general formalism of the relation between geometric phases produced by circularly 
evolving interacting spin systems and their criticality behavior is presented. This opens up 
the way for the use of geometric phases as a tool to study regions of criticality without 
having to undergo a quantum phase transition. As a concrete example a spin- 1/2 chain 
with XY interactions is presented and the corresponding geometric phases are analyzed. 
The generalization of these results to the case of an arbitrary spin system provides an 
explanation for the existence of such a relation. 
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1. Introduction 

A few conceptual advances in quantum physics have been as exciting or as broadly studied 
as geometric phases. They lie within the heart of quantum mechanics giving a surprising 
connection between the geometric properties of evolutions and their dynamics. Histori- 
cally, the first such effect was presented by Aharonov & Bohm (1959), where the quantum 
state of a charged particle acquires a phase factor when it moves along a closed path in the 
presence of a magnetic field. Since then a series of similar effects have been considered 
and experimentally verified. A more intriguing case where a quantum state is changed in a 
circular fashion was studied by Berry (1984) which led to the generation of the celebrated 
Berry phase. This effect does not need the presence of electromagnetic interactions and 
due to its abstract and general nature has found many applications (Shapere et al. 1989; 
Bohm et al. 2003). 

A characteristic that all non-trivial geometric evolutions have in common is the pres- 
ence of non-analytic points in the energy spectrum. At these points the state of the system 
is not well defined due to their degenerate nature. One could say that the generation of a 
geometric phase is a witness of such singular points. Indeed, the presence of degeneracy 
at some point is accompanied by curvature in its immediate neighborhood and a state that 
is evolved along a closed path is able to detect it. These points or regions, are of great in- 
terest to condensed matter or molecular physicists as they determine, to a large degree, the 
behavior of complex quantum systems. The geometric phases are already used in molec- 
ular physics to probe the presence of degeneracy in the electronic spectrum of complex 
molecules. Initial considerations by Herzberg and Longuet-Higgins (1963), revealed a sign 
reversal when a real Hamiltonian is continuously transported around a degenerate point. Its 
generalization to the complex case was derived by Stone (1976) and an optimization of the 
real Hamiltonian case was performed by Johansson & Sjoqvist (2004). 

Geometric phases have been associated with a variety of condensed matter and solid 
state phenomena (Thouless et al. 1982; Resta 1994; Nakamura & Todo 2002; Ryu & Hat- 
sugai 2006). Nevertheless, their connection to quantum phase transitions has only been 
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shown recently by Carollo & Pachos (2005). It was farther elaborated by Zhu (2005), 
where the critical exponents were evaluate from the scaling behavior of geometric phases, 
and by Hamma (2006), who showed that geometric phases can be used as a topological 
test to reveal quantum phase transitions. In essence, quantum phase transitions describe the 
abrupt changes on the macroscopic behavior of a system resulting from small variations of 
external parameters. These critical changes are caused by the presence of degeneracies in 
the energy spectrum and are characterized by long range quantum correlations. This is an 
exciting area of research that considers a variety of effects such as the quantum Hall effect 
and high T c superconductivity. 

Here we exploit geometric phases as a tool to probe quantum phase transitions in many- 
body systems. This provides the means to detect, not only theoretically, but also experimen- 
tally the presence of criticalities. Apart from the academic interest this approach may have 
certain advantages. In particular, the geometric evolutions do not take the system through 
a quantum phase transition. The latter is hard to physically implement as it is accompanied 
by multiple degeneracies that can take the system away from its ground state. Hence, they 
provide a way to probe criticalities in a physically appealing way. Moreover, the geometric 
phases provide a non-local object that might be useful to probe critical phenomena that 
are undetectable by local order parameters. The latter consist an exciting field of current 
research (Wen 2002). 

As an explicit example we employ the one dimensional XY model in the presence of a 
magnetic field. This model is analytically solvable and it offers enough control parameters 
to support geometric evolutions. By explicit calculations we observe that an excitation of 
the model obtains a non-trivial geometric phase if and only if it circulates a region of crit- 
icality. The generation of this phase can be traced down to the presence of degeneracy of 
the energy at the critical point in a similar way used in molecular systems. The geometric 
phase can be used to extract information on the critical exponents that completely charac- 
terize the critical behavior. A generalization of the results to the case of an arbitrary spin 
system is demonstrated. Finally, a physical implementation of the XY model and their cor- 
responding geometric evolutions is proposed with ultra-cold atoms superposed by optical 
lattices (Pachos & Rico 2005). The independence of the generated phase from the number 
of atoms, its topological nature and its resilience against control errors makes the proposal 
appealing for experimental realization. 

2. Geometric phases 

Historically, the definition of geometric phase was originally introduced by Berry in a 
context of closed, adiabatic, Schrodinger evolutions. What Berry showed in his seminal 
paper (Berry 1984) was that a quantum system subjected to a slowly varying Hamiltonian 
manifests in its phase a geometric behavior. 

Let us summarize the derivation of the Berry phase. Consider a Hamiltonian H(\) 
depending on some external parameters A = (Ai, A 2 , . . . , A m ), and suppose that these 
parameters can be varied arbitrarily inside a parameter space A4. Assume that for each 
value of A the Hamiltonian has a completely discrete spectrum of eigenvalues, 

H(X) |n(A)> = E n (X) |n(A)> , (2.1) 

where |n(A)) and E n {\) are eigenstates and eigenvalues, respectively, of H(X). Suppose 
that the values of A change slowly, along a smooth path in M. Under the adiabatic approx- 
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imation, a system initially prepared in an eigenstate |n(A)) it remains in the corresponding 
instantaneous eigenspace. 

In the simplest case of a non-degenerate eigenvalue, the evolution of the eigenstate 
is specified by the spectral decomposition (2.1) up to a phase factor. This phase factor 
can be evaluated by solving the Schrodinger equation under the constraint of the adiabatic 
approximation, yielding 

m) n ) = e' 113 exp U £ A ■ dx\ |n(A(t))) , (2.2) 

where j3 = —i f Q T E n (t)dt is the usual dynamical phase, and the extra phase factor is 
the geometric phase. This phase has the form of a path integral of a vector potential A 
(analogous to the electromagnetical vector potential) called the Berry connection, whose 
components are 

Ai = i(n(\)\-^\n(\)). (2.3) 

Berry was the first recognizing that this additional phase factor has an inherent geometrical 
meaning: it cannot be expressed as a single valued function of A, but it is a function of the 
path followed by the state during its evolution. Surprisingly, the value of this phase depends 
only on the geometry of the path, and not on the rate at which it is traversed. Hence the 
name "geometric phase". 

The simplest, but still significant example of geometric phase, is the one obtained for 
a two-level system, such as a spin- 1/2 particle in the presence of a magnetic field. Its 
Hamiltonian is given by 

H(0, 0) = B(0, <t>) ■ a = \B\U(6, cb)a z U^6, <j>), (2.4) 

where (9, 4>) determine the orientation of the magnetic field, U(8, <j)) = e - l <r z <i>/ 2 e - l <r ve / 2 
is a SU(2) transformation which rotates the operator B • <r to the z-direction and a = 
(a x , <T y ,a z ) is the vector of Pauli's operators, given by 




(2.5) 
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With this parametrization the Hamiltonian can be represented as a vector on a sphere, 
centered in the point of degeneracy of the Hamiltonian (|B| = 0), as seen in Figure 1. 

For 8 — <j> — we have that U = 1 and the two eigenstates of the system given 
by |+) = (1>0) T and |— ) = (0, 1) T with corresponding eigenvalues E± = ±|B|/2. 
Let us consider the evolution resulting when a closed path C is spanned adiabatically on 
the sphere. Following the previous general consideration it is easy to show that the Berry 
connection components corresponding, e.g. to the |+) state, are given by 

A e = 0, A = i(l-cos0) (2.6) 
that leads to the geometric phase 

<j>= <£ A-dr=\ [ sm6d6d<t>=^. (2.7) 
Jc 2 j£(0,0) 2 

The geometric phase that corresponds to the |— ) state is given by 7(C) = — f2/2.Herefi = 
J J s sin9d8d(f) is the solid angle enclosed by the loop, as seen from the degeneracy point. In 
this expression the geometric origin of the geometric phase <j> is evident. Its value depends 
only on the way in which these parameters are changed in relation with the degeneracy 
point of the Hamiltonian. 

A particularly interesting case is the one in which the Hamiltonian can be casted in a 
real form, corresponding to 9 = w/2. In this case the phase becomes (j) — tt reproducing 
the sign change of the eigenstate, when it circulates a point of degeneracy (|B| = 0), in 
agreement with Longuet-Higgins theorem. 

3. The XY model and its criticality 

In order to illustrate the connection between geometric phases and critical spin systems 
we shall consider the concrete example of a chain of spin- 1/2 particles subject to XY 
interactions. This is a one dimensional model with nearest neighbors spin-spin interactions, 
where we allow the presence of an external magnetic field oriented along the z-direction. 
The Hamiltonian is given by 

l=-M ^ ' 

where are the Pauli matrices at site I, 7 is the x-y anisotropy parameter and A is the 
strength of the magnetic field. This model was first explicitly solved by Lieb et al. (1961) 
and by Katsura (1962). Since the XY model is exactly solvable and still presents a rich 
structure it offers a benchmark to test the properties of geometric phases in the proximity 
of criticalities. 

In particular, we are interested in a generalization of Hamiltonian (3.1) obtained by 
applying to each spin a rotation with angle (f> around the ^-direction 

M 

H{4>) = U{4>)HU\<j)) with U(tf>)= Yl e'*'"*/ 2 , (3.1) 

l=-M 

in the same way as we did for the single spin- 1/2 particle. The family of Hamiltonians 
generated by varying (f> is clearly isospectral and, therefore, has the same energy spec- 
trum as the initial Hamiltonian. In addition, due to the bilinear form of the interaction term 
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we have that H{qb) is 7r-periodic in <p. The Hamiltonian H(<p) can be diagonalized by a 
standard procedure based on the Jordan-Wigner transformation and the Bogoliubov trans- 
formation (Carollo & Pachos 2005). From this procedure one can obtain the ground state, 
\g), which is given by 



\9) = LI ( cos i \°)k \0lk~ ^ sin y |l) fc |1> fc ) , (3.2) 

®fc>0 



where |0) fc and |l) fe are the vacuum and single fermionic excitation of the k-th momentum 
mode. The angle 6k is defined by cos Ok = /A/- with = cos 2 j^- — A and the energy 

gap above the ground state is given by = + 7 2 sin 2 It is remarkable that 
by inspection of (3.2) the ground state can be interpreted as the direct product of N spins 
each one having its own orientation given by the direction {2<p, Ok). As we will see in the 
next section this observation will make the evaluation of the ground state geometric phase 
a simple task. Before investigating the geometric properties of the XY model we will first 
consider the behavior of the spectrum as a function of the external parameters 7, A and <p. 

Let us first review the concept of quantum phase transitions. A many body system, 
driven by a parameter g, undergoes a quantum phase transition at a point g = g c when 
the energy density of the ground state at g = g c is non-analytic. This point is associated 
with a crossing or an avoiding of the energy eigenvalues. It is characterized by either a 
discontinuity in the first derivative of the ground state energy density (first-order phase 
transition) or by discontinuity or divergence in the second derivative of the ground state 
energy density (second-order quantum phase transition) assuming that the first derivative 
is continuous. In particular, the energy gap A between the ground and the first excited states 
vanishes like A oc \g — g c \ zl/ as g approaches g c creating a point of degeneracy that we will 
call critical. Moreover, the non-analyticity of the energy eigenvalues is related to abrupt 
changes in the ground state properties. As these transitions occur at zero temperature they 
are driven purely by quantum fluctuations (Sachdev 2001). It can be shown that the length 
of the associated quantum correlations, £, diverges like oc \g — g c \ v as g approaches the 
critical point g c . The parameters z and v are called the critical exponents and their values 
are universal, independent of most of the microscopic details of the system Hamiltonian. 

For the case of the XY model one can identify the critical points by finding the regions 
where the energy gap Ak vanishes. Indeed, there are two regions in the A, 7 space that 
are critical. When 7 = we have Ak — for — 1 < A < 1, which is a first order phase 
transition with an actual energy crossing and critical exponents z = 2 and v = 1/2. The 
other critical region is given by A = ±1 where we have A^ = for all 7. This is a second 
order quantum phase transition with energy level avoiding. When 7 = 1 and A = ±1 we 
obtain the Ising critical model with critical exponents z = 1 and v = 1. 

Finally, let us consider the criticality behavior of the rotated Hamiltonian, H{<p). The 
energy gap, A fc does not depend on the angle 0, as this parameter is related to an isospectral 
transformation. Hence, the criticality region for the rotated Hamiltonian, H(4>), is obtained 
just by a rotation of the critical points of the XY Hamiltonian around the A axis. This is 
illustrated in Figure 2, where the Ising type criticality corresponds now to two planes at 
A = 1 and A = — 1 and the XX criticality remains along the A axis for |A| < 1. 
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Figure 2. (a) The regions of criticality of the XY Hamiltonian are presented as a function of the 
parameters A and 7 and (b) the corresponding ones for the Hamiltonian H(cj)) where cj> parameterizes 
a rotation around the A axis. Possible paths for the geometric evolutions are depicted spanned by 
varying the parameter <j). 




Figure 3. The geometric phase corresponding to the ground state (a) and the relative one between 
the ground and first excited state (b) as a function of the path parameters A and 7. Each point of the 
surface corresponds to the geometrical phase for a path that is spanned by varying (f> from to 7r 
for certain A and 7. The values of the geometric phase corresponding to the loops Ci, C2 and C3 in 
Figure 2 are also indicated. 



4. Geometric phases and XY criticalities 

Figure 2 depicts the critical points of the XY model. Now we are interested in spanning 
looping trajectories in the space comprising of the Hamiltonian parameters A, 7 and <j>. 
The aim is to determine the geometric evolutions corresponding to these paths and relate 
them to the criticality regions. A special family of paths is of interest that circulates the A 
axis just by varying the parameter from zero to w. Indeed, these paths enclose the XX 
criticality only when — 1 < A < 1. As we shall see in the following it is possible to evaluate 
the corresponding geometric phases of the ground and the first excited states as a function 
of A and 7. 

Using the standard formula it is easy to show that the geometric phase of the ground 
state \g) is given by 

^ 9 = ~ l I (5 ' = 5> (1 ~ cos6k) - (4 ' 1} 
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This result can be understood by considering the form of \g), which is a tensor prod- 
uct of states, each lying in the two dimensional Hilbert space spanned by |0) fc |0)_ fe and 
|1), For each value of k(> 0), the state in each of these two-dimensional Hilbert 

spaces can be represented as a Bloch vector with coordinates (2<f), 6 k ). A change in the 
parameter <fi determines a rotation of each Bloch vector about the z-direction. A closed 
circle will, therefore, produce an overall phase given by the sum of the individual phases 
as given in (4.1) and illustrated in Figure 3(a). 

Of particular interest is the relative geometric phase between the first excited and 
ground states given by the difference of the geometric phases acquired by these two states. 
The first excited state is given by 

K> Hi) fe0 \oi kB II ( cos y l°>* \°U- ^ sin y I 1 ). I 1 )-*) - < 4 - 2 ) 

with k corresponding to the minimum value of the energy Afe. The behavior of this state is 
similar to a direct product of only N — 1 spins oriented along (2</>, 6> fe ) where the state of the 
spin corresponding to momentum ko does not contribute any more to the geometric phase. 
Thus the relative geometric phase between the ground and the excited states becomes 

(Peg = fe - <Pg = -7r(l - COS 6 ko ) (4.3) 

In the thermodynamical limit (N — > oo), (f> eg takes the form 

0, for |A| > 1 - 7 2 

where the condition |A| > 1 — 7 2 constrains the excited state to be completely oriented 
along the z-direction resulting in a zero geometric phase. As can be seen from Figure 3(b), 
the most interesting behavior of (p eg is obtained in the case of 7 small compared to A. In 
this case <p eg behaves as a step function, giving either it or phase, depending on whether 
|A| < 1 or |A| > 1, respectively. This behavior is precisely determined from whether the 
corresponding loop encloses a critical point or not and can be used as a witness of its 
presence. In particular, in the |A| < 1 — case the first term corresponds to a purely 
topological phase, while the second is a geometric contribution. Indeed, the first part gives 
rise to a phase which depends solely on the topological character of the trajectory traced by 
the (20, 9k) coordinates. In particular if n circulations are performed then the topological 
phase is nir, where n is the winding number. The second term is geometric in nature and 
it can be made arbitrarily small by tuning appropriately the couplings A or 7. This idea is 
illustrated in Figure 4, where the energy surface of ground and first excited state is depicted. 
The point of degeneracy is the intersection of the two surfaces. This is the point where 
the energy density is not analytical. Consider the case of a family of loops converging 
to a point. In the trivial case where the limiting point does not coincide to a degeneracy, 
the corresponding geometric phase converges to zero. If instead, the degeneracy point is 
included, the geometric phase tends to a finite value (Hamma 2006). 

To better understand the properties of the relative geometric phase, we focus on the 
region of parameters with 7 <C 1. In this case, it can be shown (Carollo & Pachos 2005) 
that the Hamiltonian, when restricted to its lowest energy modes, can be casted in a real 
form and, for |A| < 1, its eigenvalues present a conical intersection centered at 7 = 0. It 
is well known that when a closed path is taken within the real domain of a Hamiltonian, a 
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Figure 4. The conical intersection between the two lowest energy levels of the Hamiltonian as a 
function of its parameters. A contractible loop, i.e. a loop that can be continuously deformed to a 
point of the domain, produces a zero geometric phase. A non-trivial geometric phase is obtained 
from non-contractible loops. 



topological phase shift ir occurs only when a conical intersection is enclosed. In the present 
case, the conical intersection corresponds to a point of degeneracy where the XX criticality 
occurs and it is revealed by the topological term in the relative geometric phase tp eg . It is 
worth noticing that the presence of a conical intersection indicates that the energy gap 
scales linearly with respect to the coupling 7 when approaching the degeneracy point. This 
implies that the critical exponents of the energy, z, and of the correlation length, v, satisfy 
the relation zv = 1 which is indeed the case for the XX criticality (Sachdev 2001). In the 
following we shall see that geometric phases are sufficient to determine the exact values 
of the critical exponents and thus provide a complete characterization of the criticality 
behavior. 



We shall show here that the vacuum expectation value of a hermitian operator, O, can 
be written in terms of a geometric phase. This is a rather general result that can be used 
to study critical models, usually probed by the behavior of vacuum expectation values 
of observables, just by considering geometric phases. We assume, first, that O does not 
commute with the Hamiltonian, a requirement satisfied for the case of a non-degenerate 
spectrum and, second, that O can transform the ground state in a cyclic fashion. The latter 
provides the looping trajectories of the geometric evolutions. 

To show that let us extend the initial Hamiltonian, H , of the model in the following 
way 



5. The general case 



H(X) =H a + XO 



(5.1) 



Turning to the interaction picture with respect to O we obtain 



flim(A) = U(\t)H U\\t) 



(5.2) 
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where U(Xt) = cxp(-iXOt). From the cyclicity requirement there exists time T such 
that the unitary rotation U(XT) takes the ground state \ip) to itself, i.e. U(XT) \ip) = \ip). 
Hence, the desired cyclic evolution is obtained by a rotation generated by O. The geometric 
phase that result from the cyclic evolution is given by (2.2) and, thus, we have 



Hence, the expectation value of an operator that can generate circulations of the ground 
state is expressible with respect to a geometric phase. 

One can easily verify this relation for the simple case of a spin- 1/2 particle in a mag- 
netic field. When the direction of the magnetic field is changed adiabatically and isospec- 
trally then the state of the spin is guided in a cyclic path around the z-direction. The gen- 
erated phase is given by ip — w(l — cos 8), where 9 is the fixed direction of the magnetic 
field with respect to the z-direction. On the other hand, one can easily evaluate that the 
expectation value of the operator (1 — er z )/2 that generates the cyclic evolution is given 
by (ip\ (1 — ct z )/2 \tp) = (1 — cos#)/2 which verifies relation (5.3) as for this example 
AT = 2tt. 

This connection has far reaching consequences. It is expected that intrinsic properties of 
the state will be reflected in the properties of the geometric phases. The latter, as they result 
from a physical evolution can be obtained and measured in a conceptually straightforward 
way. Here, we are interested in employing geometric phases to probe critical phenomena 
of spin systems. Indeed, from the particular example of the XY model we saw that the 
presence of critical points can be detected by the behavior of specific geometric evolutions 
and the corresponding critical exponents can be extracted. This comes as no surprise as 
one can choose geometric phases that correspond to the classical correlations of the sys- 
tem (expectation values, e.g. of erf erf,) from where the correlation length and the critical 
behavior can be obtained. 

Let us apply this idea to the XY model we studied earlier. There the rotations are 
generated by the operator O = ^i&i- Hence, the resulting geometric phase is proportional 
to the total magnetization 



It is well known (Sachdev 2001) that the mangetization M z can serve as an order parameter, 
from which one can derive all the critical properties of the XY model just by considering its 
scaling behavior. Indeed, Zhu 2006 has considered the scaling of the ground state geometric 
phase of the XY model from where he evaluated the Ising critical exponents. As it has been 
shown here this is a general property that can be applied to any critical system. 



This construction, apart from its theoretical interest, offers a possible experimental method 
to detect critical regions without the need to cross them. When a physical system is forced 
to go through a critical region then excited states may become populated due to the vanish- 
ing energy gap, thus undermining the identification of the system state. Hence, being able 
to probe the critical properties of a physical systems just by evolving it around the critical 
area is of much interest to experimentalists as the energy gap can be kept to a finite value. 

In particular, we shall implement this model with optical lattices a system that has 
proven versatile in the field of quantum simulations. To this end, consider two bosonic 



V = XT{^\0\^) 



(5.3) 




(5.4) 



6. Physical implementation with optical lattices 
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species labelled by a = a, b that can be given by two hyperfine levels of an atom. Each 
one can be trapped by a laser field configured as a standing wave that is heavily detuned 
from any atomic transitions. Thus, the atom acts as a dipole in the presence of a periodic 
sinusoidal trapping potential that can generate one, two or three dimensional lattices. Here 
we will restrict to the case where the atoms in an arbitrary superposition of state a and b 
are confined in a one dimensional array by the help of two in-phase optical lattices. The 
tunnelling of atoms between neighboring sites is described by 

vr = -E( J - a l a c+i)-+ Rc -) t 6 - 1 ) 

la 

where a a and a\ are the annihilation and creation operators of particles o and J a their 
corresponding tunnelling coupling. When two or more atoms are present in the same site, 
they experience collisions given by 

# (0) = E -f a l4a'^ (6-2) 
laa' 

where U aa / are the collisional couplings between atom species a and a' . We shall consider 
the limit J <C U where the system is in the Mott insulator regime with one atom per lattice 
site (Kastberg et al. 1995; Raithel et al. 1998). In this regime, the effective evolution is 
obtained by adiabatic elimination of the states with a population of two or more atoms per 
site, which are energetically unfavorable. Hence, to describe the Hilbert space of interest, 
we can employ the pseudospin basis of | f) = \nf = l,n\ = 0) and | J.) = \nf — 
0, n\ = 1), for lattice site I, and the effective evolution can be expressed in terms of the 
corresponding Pauli (spin) operators. 

It is easily verified by perturbation theory that when the tunnelling coupling of both 
atomic species is activated, the following exchange interaction is realized between neigh- 
boring sites (Kuklov & Svistunov 2003; Duan et al. 2003) 

In order to create an anisotropy between the x and y spin directions, we activate a tunnelling 
by means of Raman couplings (Duan et al. 2003). Application of two standing lasers L\ 
and L2, with zeros of their intensities at the lattice sites and with phase difference <fi, can 
induce tunnelling of the state |+) = (e _ ^/ 2 |a) + e^/ 2 |6))/\/2. The resulting tunnelling 
term is given by V c — J c J2i( c l c i+i + H.c), where c; is the annihilation operator of 
|+) state particles. The tunnelling coupling, J c , is given by the potential barrier of the 
initial optical lattice superposed by the potential reduction due to the Raman transition. 
The resulting evolution is dominated by an effective Hamiltonian given, up to a readily 
compensated Zeeman term, by 

^ l> ab j 

where U{4>) was defined in (3.1). Combining the rotationally invariant Heisenberg inter- 
action H\ with H 2 gives the rotated XY Hamiltonian described by equation (3.1), where 
the parameter 7 is given by J^/(2eU a b) and e = (2J a Jb + Jc/2)/U a b is the overall en- 
ergy scale multiplying the Hamiltonian (3.1). The magnetic field term of is easily 
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produced by a homogeneous and heavily detuned laser radiation. If the radiation has am- 
plitude £1 and detuning A then the magnetic coupling is given by A = £1 2 /A. Finally, the 
angle <j> of the rotated XY Hamiltonian is given by the phase difference, cj>, of the lasers Li 
and L 2 - Hence, the complete control of the rotated XY Hamiltonian can be obtained by the 
optical lattice configuration presented here and the geometric evolutions can be performed 
by varying the phase <f) from to ir. 

7. Conclusions 

In this article, we have presented a method that theoretically, as well as experimentally, 
allows for the detection of regions of criticality through the geometric phase, without the 
need for the system to experience quantum phase transitions. The latter is experimentally 
hard to realize as the adiabaticity condition breaks down at the critical point and the state of 
the system is no longer faithfully represented by the ground state. Alternatively, the finite 
energy gap that is present along the circular procedure is sufficient to adiabatically 
prevent unwanted transitions between the ground state and the excited ones, even in the 
thermodynamical limit. 

The origin of the geometric phase can be ascribed to the existence of degeneracy points 
in the parameter space of the Hamiltonian (Hamma 2006). Hence, a criticality point can 
be detected by performing a looping trajectory around it and detecting whether or not a 
non-zero geometric phase has been generated. For the case of the XY model the topologi- 
cal nature of the resulting phase pinned to the value, f eg w n, is revealed by its resilience 
with respect to small deformations of the loop. This characteristic results from the conical 
intersection structure of the potential surfaces that is equivalent to having the critical ex- 
ponents satisfying zv = 1. In addition, the critical exponents can be evaluated by scaling 
arguments of the geometric phases (Zhu 2005). Hence, additional information about the 
critical exponents can be deduced from the topological nature and the exact value of the 
geometric phase. Moreover, topological phases are inherently resilient against control er- 
rors, a property that can be proved to be of a great advantage when considering many-body 
systems. Such a study can be theoretically performed on any system which can be analyt- 
ically elaborated such as the case of the cluster Hamiltonian (Pachos & Plenio 2004), or 
exploited numerically when analytic solutions are not known. The generalization of these 
results to a wide variety of critical phenomena and their relation to the critical exponents 
is a promising and challenging question which deserves extensive future investigation. 
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